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ABSTRACT 

A  laterally  averaged  one  dimensional  P  and  S  wave  velocity  structure  is 
obtained  for  Silent  Canyon  Caldera,  located  beneath  Pahute  Mesa  at  the  Nevada 
Test  Site.  The  velocity  models  are  derived  from  a  linearized  damped  least 
squares  travel  time  inversion  for  slowness  gradient  perturbations.  A  total  of  72  P 
wave  and  20  S  wave  travel  times  generated  from  13  different  explosions,  and 
obtained  at  scattered  locations  above  the  caldera,  are  used  in  the  inversion.  The 
structure  is  modeled  using  linear  velocity  gradients.  The  inversion  produces  a 
model  which  increases  in  velocity  from  1.50  km/sec  at  the  surface  to  1.62  km/sec 
at  2.5  km  depth  for  P  waves,  and  from  0  99  km/sec  at  the  surface  to  2.35  km/sec 
at  2.5  km  depth  for  S  waves.  The  P  and  S  wave  models  produce  travel  time  resi¬ 
duals  which  are  typically  less  than  0.1  and  0.2  seconds,  respectively,  which  is 
within  observational  error.  However,  P  wave  amplitude  ratios  and  waveform 
characteristics  of  data  collected  from  7  explosions  buried  near  625  m  depth  indi¬ 
cate  the  presence  of  lateral  velocity  variations  within  the  caldera  which  may  be 
significant  to  nuclear  source  studies. 


Introduction 


The  Silent  Canyon  Caldera,  located  within  Pahute  Mesa  in  the  north-west  corner  of  the 
Nevada  Test  Site,  has  been  and  continues  to  be  the  site  of  numerous  underground  nuclear  detona¬ 
tions.  The  ground  motion  recordings  generated  by  these  explosions  provide  the  necessary  data  for 
a  study  of  (he  physics  of  the  explosive  source  as  a  generator  of  seismic  energy.  When  inferring 
properties  of  the  seismic  source  from  observed  ground  motion,  one  is  faced  with  a  fundamental 
problem  of  seismology,  namely  that  the  net  recorded  motion  represents  a  coupling  between  the 
source  signature  and  the  response  to  the  physical  structure  through  which  the  seismic  energy  pro¬ 
pagates.  To  accurately  determine  the  characteristics  of  the  source,  the  velocity,  density,  and 
attenuation  structure  must  be  known.  The  intent  of  this  study  is  to  estimate  an  improved  shal¬ 
low  velocity  structure  of  Silent  Canyon  Caldera  in  the  hopes  that  seismic  sources  generated  there 
may  be  better  resolved  in  subsequent  studies  of  ground  motion  data. 

There  have  been  a  number  of  published  velocity  models  for  Pahute  Mesa,  e.g.,  Carrol 
(1966),  Hamilton  and  Mealy  (1969),  Helniberger  and  Hadley  (1981),  Ilartzell,  et  al.  (1983),  and 
Stump  and  Johnson  (1981).  These  models  are  based  primarily  on  a  limited  amount  of  drill  hole 
and  travel  time  data.  The  velocity  model  to  be  presented  in  this  paper  uses  a  much  larger  set  of 
travel  times  in  a  linearized  inversion  for  slowness  gradient  perturbations  from  an  initial  model. 
Additionally,  a  suite  of  near-field  waveform  data  have  been  analyzed  for  indications  of  coherent 
reflectors  within  the  caldera.  These  data  have  been  collected  at  more  or  less  random  locations 
above  the  caldera  (Figures  la  and  lb).  Reflection  or  refraction  surveys  specifically  aimed  at  eluci¬ 
dating  the  two  or  three-dimensional  structure  of  the  caldera  have  not  yet  been  attempted.  There¬ 
fore,  in  the  absence  of  a  more  extensive  data  set,  we  will  offer  a  one  dimensional  velocity  structure 
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which  averages  through  any  lateral  variations.  Previous  velocity  models  have  likewise  made  this 


simplification.  In  an  attempt  to  assess  the  validity  of  this  approximation,  we  will  qualitatively 


examine  the  possible  effect  of  lateral  heterogeneity  on  the  suite  of  waveform  data.  Following  this 


will  be  the  results  of  the  one  dimensional  travel  time  inversion. 


Geologic  Setting  and  Lateral  Heterogeneity 


The  Silent  Canyon  Caldera  is  located  within  a  pyroclastic  volcanic  center  located  in  the 


eastern  portion  of  Pahute  Mesa.  The  center  was  volcanically  active  for  at  least  1  m.y.  in  the  late 


Miocene  (13.1  -  1-1.8  m.y.)  during  which  time  it  collapsed  and  became  partially  filled  with  self- 


derived  tuffs  and  lava  flows  (Noble,  et.  al.  1968).  Subsequent  to  the  collapse  and  general  volcanic 


activity,  the  caldera  became  almost  completely  buried  with  volcanics  erupted  from  younger 


nearby  centers.  The  surface  outline  of  the  caldera  shown  in  Figures  la  and  lb  is  for  the  most 


part  inferred  from  gravity,  surface  geological  mapping,  and  subsurface  data. 


The  availability  of  limited  borehole  data  throughout  the  caldera  has  made  it  possible  for 


geologists  to  construct  stratigraphic  profiles.  Volcanic  sediments  have  been  locally  mapped  down 


to  -1  km  depth.  The  overall  interpretation  by  Orkild,  et.  al.  (1968)  indicates  three  basic  types  of 


structural  lateral  heterogeneities.  These  are  (l)  an  abundance  of  basin  and  range  normal  faulting, 


(2)  interfingering  of  rock  units,  and  (3)  non-uniform  collapse  of  the  caldera  itself.  These  structural 


variations  are  reflected  in  down  hole  density  measurements  which  vary  from  about  1.9  to  2.5 


gm/cc  at  a  given  depth  down  to  at  least  1.5  km  among  boreholes  throughout  the  caldera  (Ifealy, 


.V  .V.VLV.s 
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Our  one-dimensional  model  cannot  deal  with  such  irregularities.  However,  in  seismic  source 
studies  it  is  important  to  at  least  qualitatively  understand  the  effect  of  lateral  heterogeneities  on 
waveform  data.  To  this  end,  we  now  consider  the  near  field  ground  motion  data  from  seven 
explosions,  each  buried  near  625  m  depth,  and  examine  the  effect  that  lateral  variations  may  have 
on  these  data.  Combining  explosions  with  similar  source  depths  allows  us  to  interpret  the 
waveforms  as  if  (hey  were  all  generated  from  a  single  source.  The  explosions  used  here  are  Chan¬ 
cellor,  Colwick,  Emmenthal,  Farm,  Ilarzer,  Pipkin,  anti  Serpa  (see  Figure  la.  and  Table  la) 
Most  stations  recorded  all  three  components  of  ground  acceleration,  yielding  a  total  of  4-1  record¬ 
ings  per  component.  The  recordings  from  these  scattered  stations  are  contracted  onto  a  one¬ 
dimensional  line  and  plotted  in  Figures  2  to  4  as  a  function  of  distance  away  from  their  respective 
explosion.  Each  trace  is  normalized  to  its  own  maximum  amplitude.  For  reasons  to  be  discussed 
later,  the  traces  have  been  shifted  slightly  in  time. 

The  actual  first  arrival  times  of  these  waveforms,  as  measured  off  the  vertical  component, 
are  plotted  in  Figure  5.  A  fair  amount  of  scatter  is  evident,  ranging  up  to  0.3  seconds.  If  the  cal¬ 
dera  were  laterally  homogeneous,  then  the  scatter  of  points  could  only  be  due  to  timing  inaccura¬ 
cies.  We  estimate  the  instrumental  timing  errors  of  these  data  to  be  typically  no  greater  than 
about  0.05  seconds,  though  the  scatter  in  the  Colwick  times  suggests  possible  timing  errors  up  to 
a  tenth  of  a  second.  The  picking  error  is  typically  about  0.01  seconds  for  these  P  arrivals.  The 
scatter  may  also  in  part  be  due  to  differences  in  station  elevations.  The  average  height  of  the  sta¬ 
tions  above  the  elevation  of  the  lowest  station  is  132  m,  the  maximum  being  225  m.  For  a  near 
surface  velocity  of  1.5  km,  this  would  result  in  relative  delay  times  of  a  tenth  of  a  second  or  so. 
However,  an  attempt  to  correct  the  travel  times  to  a  common  datum  did  not  reduce  the  scatter 
by  any  significant  amount.  The  scatter  of  these  travel  times  beyond  the  total  timing  and  picking 


error  of  0.11  seconds  thus  suggests  that  there  exist  lateral  velocity  variations  capable  of  generat¬ 
ing  travel  time  anomalies  up  to  a  few  tenths  of  second.  The  ray  bottoming  depths  for  these  scat¬ 
tered  arrivals,  as  determined  from  our  final  model,  suggest  that  such  heterogeneities  are  present 
down  to  a  depth  of  at  least  2.5  km,  consistent  with  the  geologic  interpretation. 

An  additional  indication  of  lateral  heterogeneity  is  the  variation  in  the  ratio  of  P  wave  vert¬ 
ical  amplitude,  .-I;-,  to  P  wave  radial  amplitude,  AR  .  If  the  media  were  laterally  homogeneous 
and  each  of  the  seven  explosions  generated  similar  source  pulses,  one  would  expect  similar  varia¬ 
tions  of  Ay/AR  as  a  function  of  epicentral  distance.  The  P  wave  acceleration  amplitude  ratios 
for  these  data,  as  measured  off  the  first  quarter  cycle,  are  given  in  Figure  6.  Significantly  different 
patterns  of  Ay  / AR  between  explosions  are  apparent.  This  is  an  indication  that  waves  travelling 
the  same  epicentral  distance,  but  in  different  parts  of  the  caldera,  may  be  sampling  substantially 
different  velocities  along  their  paths. 

Some  degree  of  this  disparity  between  explosions  may  also  be  due  to  topographic  focussing 
effects.  Additionally,  many  of  the  stations  recorded  transverse  motion  arriving  coincident  with 
the  first  half  cycle  of  the  P  arrival  on  the  two  isotropic  components,  i.e.,  the  radial  and  vertical 
components  (see  Figure  -1).  The  amplitudes  of  these  arrivals  are  often  a  few  tens  of  percent  of  the 
vertical  P  wave  amplitude  and  even  as  large  as  95  %  at  one  of  the  Farm  stations.  These 
transverse  arrivals  may  be  indicative  of  lateral  diffraction  and/or  scattering  of  energy,  both  of 
which  may  affect  the  P  wave  amplitudes  on  the  vertical  and  radial  components  of  ground  motion. 
Though  the  implications  of  the  amplitude  data  are  not  entirely  conclusive,  some  observable 
degree  of  lateral  inhomogeneity  in  the  structure  is  indicated. 


Our  final  comments  on  the  subject  of  lateral  heterogeneity  within  the  caldera  are  directed  at 
the  waveforms  themselves.  By  combining  t he  traces  into  record  sections,  as  in  Figures  2  to  1,  our 
intent  was  to  search  the  suite  of  seismograms  for  indications  of  subsurface  structure,  i.e., 
reflections  passing  coherently  through  the  sections.  The  problem  of  static  time  corrections  was 
best  solved  by  fitting  a  smooth  curve  through  the  scattered  P  travel  times  in  Figure  5  and  shifting 
the  traces  on  all  three  components  by  the  amount,  necessary  to  align  the  vertical  P  wave  onsets 
with  the  smooth  curve.  The  curve  was  initially  generated  in  a  trial  and  error  fashion  by  numeri¬ 
cally  calculating  the  P  wave  travel  times  for  various  preliminary  velocity  models.  The  model 
which  appeared  to  best  average  through  the  data  also  served  as  the  initial  velocity  model  for  the 
linearized  inversion.  In  an  attempt  to  align  later  arrivals  in  the  sections  as  best  as  possible,  the 
first  arrivals  were  shifted  to  coincide  with  the  predicted  P  wave  travel  times  from  our  final  model, 
as  shown  in  Figure  5. 

Though  reflection  hyperbolas  are  not  clearly  evident,  the  sections  do  show  some  coherency 
between  the  traces  recorded  at  the  closely  packed  Colwick  and  Serpa  stations.  In  an  attempt  to 
extract  some  amount  of  structural  detail  from  such  coherency,  the  sections  were  subjected  to 
some  of  the  standard  processing  techniques  that  are  applied  to  controlled  source  reflection  data. 
For  example,  the  traces  were  filtered  at.  various  bandpasses  and  examined  for  reflections  using  a 
velocity  spectra  technique,  (Taner  and  Koehler,  1069).  However,  no  definitive  reflecting  horizons 
were  uncovered  in  this  process.  This  may  in  part  be  due  to  negligible  acoustic  impedance  con¬ 
trasts  across  stratigraphic  units.  However,  the  overall  non-coherency  between  events  and  within 
non-array  recorded  events  suggests  that  the  waves  have  traveled  through  and  have  been  reflected 
off  different  lateral  heterogeneities  within  the  caldera.  The  presence  of  near  surface  hetero¬ 
geneities  are  suggested  by  the  non-uniformity  of  surface  waves  across  the  sections. 
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These  waveform  data,  as  well  as  P  wave  travel  times  and  amplitude  measurements,  are 
fairly  strong  indicators  that  lateral  variations  of  structure  within  Silent  Canyon  Caldera  may  be 
significant.  However,  without  good  quality  reflection  or  refraction  data,  one  must  settle  for  an 
averaged  one  dimensional  velocity  model. 


Travel  Time  Inversion 


In  an  attempt  to  make  more  quantitative  use  of  the  available  travel  times,  we  have  imple¬ 
mented  a  least  squares  travel  time  inversion.  Given  an  initial  velocity  model  reasonably  close  to 
the  true  model,  one  can  make  use  of  travel  time  data  to  refine  the  model  parameters,  bringing 
them  closer  to  a  model  which  better  fits  the  travel  times.  Work  by  Backus  and  Gilbert  (1969) 
and  Pavlis  and  Booker  (1980)  describe  the  method  of  velocity  model  refinement  in  detail  and 
apply  it  to  earthquake  travel  time  data.  In  our  study,  using  only  explosion  data,  we  have  the 
advantage  that  locations  and  origin  times  are  known  and  need  not  be  solved  for  in  the  inversion. 
The  problem  is  then  Amplified  and  is  outlined  below. 

For  a  group  of  non-dispersive  body  waves,  the  travel  time  from  a  hypocenter  //  to  a 
receiver  A,-  is 

T{  =  J  vds  (1) 

H^t i, 

where  «  is  the  reciprocal  velocity,  or  slowness,  and  da  is  the  element  of  path  length.  Ignoring 
second  order  variations  in  T  due  to  the  assumed  small  change  in  integration  path,  a  small  change 
in  the  slowness  6u  [z )  results  in  a  change  in  travel  time  of 

STt  =  /  6i i(z)rf»  (2) 

H^, 
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For  the  flat  layered  geometry  used  here  we  have  for  the  i‘*  travel  time  observation  with  ray 
parameter  p,  , 


6T,  =  /  6u(z)9i(z)dz 

L 


(3) 


where  J,  (;)  depends  on  the  source  depth  z//  as  follows: 


?.  (-  )  "  I 


Mil 


(M(r)-p,2)'2 

2«  (-’  ) 

(«‘(M  -  P.  T 

0 


0  <  i  <  z„ 

Zfi  <z  <  depth  of  ray  bottom  (u  =  p, 
2  >  depth  of  ray  bottom 


(4) 


/,  is  chosen  in  (3)  as  some  depth  below  the  deepest  ray.  Notice  that  the  kernel  g,(z  )  is  singular 
at  t he  turning  point  where  u  —  p,  .  This  singularity  can  be  removed  through  integrating  by 
parts,  yielding 


o 

<5T,  -  Su  (0)(7,  (0)  -  fbu' (z)0,(z)Jz 

L 


(5) 


where 


bu'  (;  )  —  -j-Su  (z  )  and  (,\(z  )  /s,(s)'M 

/ 


Discretizing  the  media  into  plane  layers  allows  one  to  solve  (5)  for  slowness  gradient  pertur¬ 
bations  b u  (;)  and  the  surface  slowness  perturbation  bu  (0)  using  standard  linear  inversion  tech¬ 
niques.  The  gradients  can  then  be  integrated  to  yield  a  revised  slowness,  or  equivalently,  velocity 
model.  The  process  can  be  continued  until  the  root  mean  square  (RMS)  travel  time  residual  is 
approximately  equal  to  the  standard  error  of  the  data.  Practically  speaking,  the  discretization  of 
the  media  will  remove  the  singularity  in  (5),  and  one  could  also  solve  for  absolute  slowness 
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perturbations  with  relative  ease.  The  choice  here  of  solving  for  the  gradient  is  twofold.  First,  it 
is  the  change  in  velocity  with  depth  rather  than  absolute  values  of  velocity  that  governs  the 
focussing  of  seismic  energy  and  thus  largely  determines  waveform  characteristics.  Secondly,  the 
integration  of  the  kernels  j,  (j)  results  in  smoother  solutions,  which  is  important  for  stable  itera¬ 
tions. 

The  above  inversion  for  slowness  gradient  perturbations  was  carried  out  using  a  damped 
singular  value  decomposition  algorithm  (O’Connell,  1086).  The  damping  provides  additional  sta¬ 
bility  in  this  typically  unstable  problem  of  velocity  inversion.  Undamped  error  and  resolution  esti¬ 
mates  were  made  using  C  =  and  R  —  VVT  ,  respectively,  where  C  is  the  covariance 

matrix,  R  is  the  resolution  matrix,  and  G=UAVT  is  the  standard  singular  value  decomposition  of 
the  kernel  matrix,  a  here  represents  the  smallest  apriori  variance  in  the  data  set,  i.e.,  the  instru¬ 
mental  timing  errors  alone.  Because  of  the  various  uncertainties  in  picking  arrival  times,  the 
equations  in  (5)  were  each  weighted  by  their  own  uncertainties  relative  to  a.  Inherent  in  this 
analysis  is  the  tradeoff  between  error  and  resolution  as  a  function  of  the  number  of  singular  values 
retained  in  the  singular  value  decomposition.  To  achieve  reasonable  error,  the  smaller  singular 
values  must  be  disposed  of  at  the  price  of  poorer  resolution.  It  is  important  to  remember  (hat  C 
and  R  are  only  linear  estimates  for  this  non-linear  problem.  The  pitfalls  of  ignoring  the  non¬ 
linear  effect  are  discussed  by  Pavlis  and  Booker  (1983). 

The  travel  time  data  which  were  used  in  the  inversion  are  listed  in  Tables  la  and  lb,  which 
correspond  to  the  source  and  receiver  locations  shown  in  Figures  la  and  lb.  Three-component 
acceleration  records  were  available  for  the  events  in  Table  la  and  the  arrival  times  of  P  and  S 
waves  were  measured  off  these  records.  S  waves,  assumed  to  have  been  generated  at  the  site  of 
the  explosion,  were  picked  off  the  transverse  component  of  ground  motion  only  when  they 
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appeared  particularly  evident.  Integrated  accelerations  aided  in  the  selection.  We  estimate  the 
overall  accuracy  of  these  S  times  to  be  approximatly  0.2  seconds.  These  data  were  supplemented 
with  the  P  travel  times  listed  in  Table  lb,  which  were  taken  from  the  studies  of  Perret  (1976)  and 
Barker,  et.  al.  (1985).  Note  that  the  events  in  Table  lb  greatly  extend  the  depth  range  of  sources 
used  in  the  inversion,  which  helps  constrain  the  velocity.  Also  note  that  the  depth  to  the  water 
table  for  the  events  listed  in  Table  la  and  lb  ranges  between  580  km  and  700  km,  with  an  aver¬ 
age  of  650  m  (Howard  and  Richardson,  1981).  Thus,  sources  both  above  and  below  the  water 
table  have  been  included.  To  compensate  for  variations  in  topography,  the  shot  depth  used  for  a 
given  shot  was  that  relative  to  the  average  elevation  of  the  stations  recording  the  event.  The 
actual  depths  of  burial  are  given  in  Tables  la  and  lb.  The  total  number  of  data  used  in  the 
inversion  was  72  P  wave  times  and  20  S  wave  times. 

The  inversion  method  described  above  requires  that  a  reasonably  good  initial  estimate  of  the 
velocity  model  be  available.  The  initial  P  wave  model  for  the  upper  700  m  was  obtained  by 
averaging  acoustic  well  log  data  derived  from  the  shot  reports  of  the  events  in  Table  la  and  from 
additional  well  data  from  Howard  and  Richardson  (198-1).  The  initial  P  wave  model  below  700  m 
was  obtained  from  a  average  of  previous  models,  making  adjustments  to  reasonably  match 
observed  travel  times  as  described  above.  Given  the  initial  P  model,  the  initial  S  wave  model  was 
obtained  from  the  work  by  Carrol  (1968),  which  relates  shear  velocity  to  compressional  velocity 
through  measurements  of  Young’s  modulus  for  Pahute  Mesa  and  nearby  Northern  Yucca  Flat  vol- 
canics.  These  initial  models  are  shown  in  Figure  7. 

For  the  purposes  of  the  inversion,  the  media  was  discretized  into  100  m  thick  layers  down  to 
a  depth  of  -1  km.  Pnlike  previous  models  which  have  characterized  Silent  Canyon  Caldera  by 
constant  velocity  layers,  we  have  used  a  linear  gradient  velocity  function.  This  should  be  more 


consistent  with  the  expected  gradual  increase  in  velocity  due  to  sediment  compaction.  Such 
behavior  is  observed  in  the  shallow  well  log  data. 

The  results  of  the  velocity  inversion  are  shown  in  Figure  7  and  Table  2.  The  inversions  for 
the  P  and  S  models  were  done  independently.  It  should  be  noted  that  the  sharp  features  of  the 
models  are  not  resolvable  by  the  travel  time  data.  Rather,  they  are  artifacts  of  the  initial  models. 
The  P  model  converged  in  3  iterations  with  an  RMS  residual  of  0.096  seconds.  The  S  model  con¬ 
verged  in  only  2  iterations  with  an  RMS  residual  of  0.112  seconds.  Both  solutions  were  obtained 
with  damping  values  equal  to  approximatley  2?c  of  the  maximum  singular  value  used  in  each 
inversion.  The  travel  time  residuals  for  these  models  are  given  in  Table  la  and  lb.  Additionally, 
Figure  5  shows  the  predicted  P  and  S  travel  time  curves  for  a  shot  buried  625  m  averaging 
through  the  data  from  shots  near  that  depth. 

The  initial  P  and  S  models  were  each  revised  down  to  2.5  km.  In  the  final  models  no  rays 
in  our  data  set  pene  .rate  below  this  depth  and  therefore  the  deeper  velocity  gradients  remain 
unchanged  from  the  initial  model.  The  total  number  of  singular  values  in  the  error  and  resolution 
analysis  is  thus  26,  including  one  for  the  surface  slowness.  Of  the  26,  only  6  could  be  retained  to 
achieve  an  acceptable  amount  of  error  for  the  P  model.  For  the  S  model  only  3  singular  values 
could  be  kept.  This  small  number  is  not  surprising  considering  the  larger  errors  of  the  S  wave 
travel  times.  Figure  8a  and  8b  show  the  magnitude  of  the  slowness  gradient  and  its  standard 
error  for  both  models.  The  statistical  fit  is  particularly  good  above  600  m.  This  is  less  than  the 
typical  depth  of  most  of  the  shots  used  in  this  study,  so  there  is  an  abundance  of  rays  sampling 
above  this  depth.  However,  the  lack  of  turning  rays  above  600  m  results  in  poor  vertical  resolu¬ 
tion  here.  This  is  shown  rather  well  in  Figures  9a  and  9b  where  elements  of  the  resolution  matrix, 
R,  for  P  and  S  wave  slowness  gradients  are  shaded  according  to  their  value.  The  poor  vertical 
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resolution  for  the  shallow  P  model  is  compensated  by  available  well  log  data.  Below  600  m  the 
increased  resolution  of  the  P  model  trades  off  with  the  increased  relative  error.  However,  due  to 
the  lack  of  deeper  shots,  a  smaller  number  of  travel  times,  and  larger  reading  errors,  much  of  the 
S  wave  model  remains  poorly  resolved. 


Summary  and  Conclusions 


The  intent  of  this  study  has  been  to  estimate  an  averaged  one-dimensional  velocity  model 
for  Silent  Canyon  Caldera.  Our  primary  data  set  has  been  acoustic  well  log  data  and  a  relatively 
large  set  of  travel  times.  These  data  have  produced  a  one-dimensional  averaged  velocity  model 
which  fits  the  observed  travel  time  data  to  within  observational  error.  For  comparison,  our  final 

w  L 

model  is  shown  in  Figure  10  with  four  recently  published  models.  A  major  difference  between  this 
and  previous  studies  is  the  introduction  of  a  continuous  velocity  model.  The  P  wave  velocity 
model  of  this  study  is  slightly  greater  than  previous  models  between  depths  of  about  0.7  and  2 
km,  but  this  is  compensated  by  the  lower  velocities  at  depths  less  than  0.3  km.  The  S  wave  velo¬ 
city  model  is  greater  than  the  results  of  most  other  studies.  Using  our  inversion  approach,  a  one¬ 
dimensional  model  could  be  even  better  constrained  with  a  broader  depth  distribution  of  sottrces 
and  more  recorded  travel  times,  particularly  those  of  S  waves.  Though  a  one  dimensional  model 
may  satisfy  travel  time  data  to  within  acceptable  error,  amplitude,  waveform,  and  geologic  data 
indicate  that  lateral  heterogeneity  within  Silent  Canyon  Caldera  is  significant.  Given  the  desire 
for  unequivocal  descriptions  of  nuclear  sources,  a  more  detailed  investigation  of  the  three  dimen¬ 
sional  seismic  structure  of  the  caldera  is  certainly  warranted  by  the  above  discussions. 
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TABLE  la 


TRAVEL  TIMES  AND  RESIDUALS 


i 


EVENT 

DEPTH 

STA. 

A 

P  TIME 

A  TP 

S  TIME 

AT* 

(m) 

(km) 

(sec) 

(sec) 

(sec) 

(sec) 

CHANCELLOR 

625 

10 

1.85 

0.790 

0  07 

1.25 

0  01 

3 

2.11 

0.786 

-0  01 

1.40 

0  03 

o 

2.67 

1.053 

0  10 

1.60 

-0  05 

0 

3.09 

1.027 

-0.04 

- 

- 

5 

4.48 

1.328 

-0  13 

2.30 

-0  27 

4 

5.52 

1.631 

-0.11 

- 

- 

9 

5.60 

1.812 

0  05 

- 

- 

6 

6  03 

1.696 

-0  17 

- 

- 

1 

9.35 

2.655 

0  00 

- 

- 

7 

9.82 

2.540 

-0.22 

- 

- 

8 

10.82 

3.103 

0.12 

5.67 

0  04 

COLYVICK 

633 

4 

4.94 

1.617 

0  02 

2.90 

1 

4.94 

1.627 

0  03 

2.92 

0.10 

o 

4.94 

1.662 

0  07 

3.00 

l  . 

7 

1.677 

0  06 

3  04 

0.17 

8 

5.03 

1.841 

0  22 

3.02 

0.15 

6 

5.1 1 

1.712 

0  07 

3.08 

0.17 

9 

5.11 

1.717 

0.08 

3.05 

0  14 

3 

5.29 

1.737 

005 

3.15 

0.15 

EMMENTILAL 

576 

4 

0.51 

0.289 

-0.04 

0  j2 

0.00 

6 

0  76 

0.354 

-0  05 

0.35 

0.02 

3 

0.89 

0.424 

-0  02 

0.42 

0.04 

FARM 

689 

4 

4  22 

1.332 

-0  04 

- 

- 

6 

4.72 

1.252 

-0.26 

- 

- 

5 

6.38 

1.796 

-0.16 

- 

- 

2 

9.73 

2.605 

-0.18 

- 

- 

HARZER 

637 

4 

3.43 

1.15 

-0  02 

- 

- 

2 

3.52 

1.03 

-0  16 

- 

- 

i 

4.65 

1.39 

-Oil 

- 

- 

6 

5.61 

1.85 

0  09 

- 

- 

8 

5.78 

1.87 

0.06 

- 

- 

5 

6.62 

2.00 

-0.03 

- 

- 

PIPKIN 

624 

1 

2.71 

0.10 

- 

- 

2 

5.72 

-■T'M* 

0.28 

- 

- 

3 

7.69 

2.56 

0.28 

- 

- 

SERPA  (+) 

630 

9 

8.63 

- 

- 

4  65 

-0  01 

7 

8.79 

- 

- 

4.70 

-0  04 

6 

8.87 

- 

- 

4.75 

-0  02 

5 

8.95 

" 

- 

4.80 

-001 

+  Gross  timing  errors.  S  times  picked  after  P  times  set  to  predicted  times. 
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TABLE  lb 


TRAVEL  TIMES  AND  RESIDUALS 


EVENT 

DEPTH 

STA. 

A 

P  TIME 

A  TP 

(m) 

(Km) 

(sec) 

(sec) 

BOXCAR  (*) 

1166 

1 

0.300 

0.47 

0.04 

3 

0.91 

0.56 

0.05 

8 

2.40 

0.90 

0.03 

12 

3.81 

1.26 

0.01 

16 

4.87 

1.51 

-0.01 

21 

7.27 

2.10 

-0.01 

31 

10.37 

2.92 

0.13 

GREEK Y  (*) 

1216 

l 

0.015 

0.430 

0.00 

o 

0.61 

0.480 

0.01 

3 

l  22 

0.600 

0.01 

HALFBEAK  (*) 

819 

1 

0.015 

0.330 

0.01 

•1 

0.46 

0.410 

0.05 

3 

0.91 

0.490 

0.03 

1 

1.37 

0.640 

0.06 

7 

2  13 

0.740 

-0.05 

5 

2  13 

0.820 

0  03 

8 

2.28 

0  880 

0  04 

6 

3  05 

1  040 

-0.01 

INLET  (#) 

819 

E 

0  01 

0  31 

-0  01 

W 

0  01 

0  31 

-0  01 

3 

0  41 

0.36 

0.01 

4 

0.82 

0  44 

0.00 

5 

1.63 

0.59 

-0  06 

6 

3.27 

1.03 

-0  08 

7 

6.53 

1.99 

0  02 

MAST  (#) 

91 1 

Z 

0  01 

0.36 

0  02 

2 

0.46 

0  40 

0.02 

3 

0.91 

0.51 

0  04 

4 

1.82 

0.71 

0.00 

5 

3  65 

1.15 

-0.06 

6 

5.47 

1.56 

-0.13 

7 

7.30 

2.03 

-0.11 

SCOTCH  (*) 

977 

1 

0  02 

0  41 

0  05 

2 

1  21 

0  64 

0  08 

3 

2.29 

0.98 

0.14 

3  A 

4.13 

1.38 

0  04 

4 

6.06 

1.89 

0.05 

*  Times  from  Perret  (1976). 

#  Times  from  Barker,  et.  al.  (1985). 
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TABLE  2 

P  AND  S  VELOCITY 


DEPTH  P  VELOCITY  S  VELOCITY 

(km)  (km/sec)  (km/sec) 
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Figure  Captions 


Figure  la.  Explosion  sites  (open  circles)  and  respective  recording  sites  (filled  circles)  for 
shots  near  625  m  depth  at  Pahute  Mesa  of  the  Nevada  Test  Site.  See  Table  la. 


Figure  lb.  Additional  explosion  sites  (open  circles)  and  respective  recording  sites  (filled  cir¬ 
cles).  See  Table  lb. 

Figure  2.  Vertical  accelerations  from  shots  shown  in  Figure  la.  Shot  depths  are  near  625  m. 
Figure  3.  Radial  accelerations  from  shots  shown  in  Figure  la.  Shot  depths  are  near  625  m. 

Figure  4.  Transverse  accelerations  from  shots  shown  in  Figure  la.  Shot  depths  are  near  625 

m. 

Figure  5.  Measured  P  wave  and  S  wave  travel  times  for  shots  buried  near  625  m.  Predicted 
times  for  a  shot  buried  625  m,  as  determined  from  our  final  model,  are  given  by  the  solid 
curves. 

Figure  6.  Ratio  of  vertical  to  radial  P  wave  acceleration  for  shots  buried  near  625  m.  Ampli¬ 
tudes  were  measured  off  the  first  quarter  cycle. 

Figure  7.  Initial  and  revised  velocity  models.  No  constraints  are  placed  on  the  models  below 
2.5  km,  thus  the  deeper  velocity  gradients  remain  unchanged  from  the  initial  model. 

Figure  8a.  Final  slowness  gradient  and  standard  error  (diagonal  elements  of  the  covariance 
matrix  C  )  for  the  P  wave  model.  The  slowness  gradient  is  actually  negative  ;  it’s  absolute 
value  is  plotted. 

Figure  8b.  Same  as  Fig.  8a.,  but  for  S  wave  model. 

Figure  9a.  Elements  of  the  resolving  matrix,  R  ,  for  the  P  wave  model.  The  values  are  are 
divided  into  six  levels  and  shaded  with  a  corresponding  density.  The  upgoing  diagonals  sig¬ 
nify  negative  values,  indicative  of  side  lobes.  The  first  row  and  column  pertain  to  the  value 
of  surface  slowness. 


Figure  9b.  Same  as  Fig.  9a.,  but  for  S  wave  model. 

Figure  10.  Final  P  and  S  wave  velocity  models  plotted  with  recently  published  models. 
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Figure  lb.  Additional  explosion  sites  (open  circles)  and  respective  recording  sites  (filled  cir¬ 
cles).  See  Table  lb. 


were  measured  off  the  first  quarter  cycle. 
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Figure  9a.  Elements  of  the  resolving  matrix,  R  ,  for  the  P  wave  model.  The  values  are  are 
divided  into  six  levels  and  shaded  with  a  corresponding  density.  The  upgoing  diagonals  sig¬ 
nify  negative  values,  indicative  of  side  lobes.  The  first  row  and  column  pertain  to  the  value 
of  surface  slowness. 
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